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Abstract 

Realistic boundaries of finite heat conductivity for thermoconvection in a 
Rayleigh-Benard setup with magnetized ferrofluids are investigated. A linear 
stability analysis of the conductive state is performed with a shooting method. 
It shows that the critical wave number is for any magnetic field stronly in- 
fluenced by the conductivity of the boundaries. Linear as well as nonlinear 
coefficients of a Ginzburg Landau amplitude equation for convection shortly 
above the onset are evaluated as functions of the magnetic Rayleigh number, 

the boundary conductivities, and the fluid Prandtl number. 
PACS: 83.80.Gv Electro- and magnetorheological fluids; 47. 20. Bp Buoyancy- 
driven instability; 47.20.Ky Nonlinearity (including bifurcation theory); 
47.65.+a Magnetohydro dynamics and electrohydrodynamics 
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I. INTRODUCTION 



Convection of magnetic fluids in the classical Rayleigh- Benard setup in the presence 
of a vertical magnetic field has been studied theoretically by Finlayson M. However, the 
experiments by Schwab et al. |2||| showed that the theoretical critical wave numbers of 
the flow structure were larger than those that were measured in the experiments done with 
horizontal plates of finite heat conductivity. In this paper we investigate how with decreasing 
conductivity of the plates linear and weakly nonlinear convection properties change. 

We consider a horizontal fluid layer of height d between horizontal plates each of thickness 
dp subjected to a homogeneous vertical (z-direction) gravitational acceleration g and a 
homogeneous vertical magnetic field. The ratio of the termal conductivities of the plates, 
Ap, and of the fluid, A F , is ( = and 5 = ^ is the plate thickness in units of half the height 
of the fluid layer. Then the vertical temperature profile in the conductive, quiescent state 
without convection is piecewise linear as shown in Fig. [I] for the case of 5 = 1. A temperature 
difference AT imposed between the outer surfaces of the plates at z = ±(| + dp) implies a 
temperature difference 

AT 

^ - TTJfc ™ 

across the fluid layer in the conductive state. We use this conductive ATp to define the 
Rayleigh number 

Ra = a9 d3ATF (1.2) 

KV 

and the magnetic Rayleigh number 

N= ^p^Ef (L3) 

K7](l +X) 

as control parameters. A negative value of AT or of ATp stands for heating of the upper 
plate and a positive one for lower plate heating. Here a is the thermal expansion coefficent, 
k the thermal diffusivity, v the kinematic viscosity, and rj the effective viscosity of the fluid. 
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Furthermore, //o is the magnetic field constant, x the magnetic susceptibility, and K is the 
pyromagnetic coefficient, which is almost proportional to the external magnetic field. 

For convenience Ra and iV are normalized by the critical Rayleigh number Ra® for onset 
of convection in a fluid bounded by plates of conductivity Ap and thickness dp in the absence 
of a magnetic field. So we introduce 

_ Ra _ . , 

II. EQUATIONS AND BOUNDARY CONDITIONS 

Henceforth we scale positions by the height of the fluid layer d, time by — , velocity by 
77, temperature by ATp, and the vertical magnetic field by d , f^ 3 . Using this scaling we get 
from the basic hydrodynamic field equations |IJ the following system of partial differential 
equations 

ifi t V 2 u> = V 4 u> + (fig + d 2 ) [{Ra + N)9 -%]-^VxVx [(v • V)v], 

+N [d y (d y 6d z <p - d z 6d y if) + d x (d x 6d z <p - d z 6d xV )} (2.1a) 
d t 9 = V 2 9 + w-{y- V)9 (2.1b) 
V 2 <p = d 2 z 6 (2.1c) 

for the deviations from the quiescent state of heat conduction in the fluid. Here w = v • e z 
is the vertical velocity field, 6 = T — T cond , <p — (H — H cond ) ■ e z , and a = - is the Prandtl 
number. The pressure was eliminated by applying twice the rotation to the momentum 
balance equation. 

The equations are solved subject to the following realistic boundary conditions: The 
plates are rigid, thus enforcing the no-slip condition 

w = d z w = at z — ±- (2.2a) 

on the velocity field. The boundary conditions on the deviations of the temperature in the 
plates, 9, and in the fluid, 9, from the conductive profile of Fig. [I] are 



d z e = (d z 6, 6 = 6 at z = ± 2 ( 2 - 2b ) 

to ensure continuity of the heat flux and of the temperature. With AT being imposed 
externally one has 

= at* = ±-(l + «f). (2.2c) 

The vertical magnetic field deviation from the magnetic field in the conductive state, tp = 
(H — H con d) ■ e z , has to fulfill the condition 

(1 + X)d g (p ± k(p = at z = ± \- ( 2 - 2d ) 

Here k is the wave number of the convective roll structure that bifurcates out of the con- 
ductive state. 

III. LINEAR ANALYSIS 

We have performed a standard linear stability analysis of the conductive state for plane 
wave perturbation fields obeying the above boundary conditions fl2.2|). To that end the 
linearized version of the field equations of the ferrofluid ( |2.1p coupled to the heat diffusion 
equation in the plates were solved by a shooting method ||. The resulting critical eigen- 
functions are plotted in Fig. ^ for different magnetic fields at 5 = 1 and ( = 5. The latter 
being the conductivity ratio realized in the experimental setup of Schwab et al. |2|||. The 
maximal amplitude of the vertical velocity w(z) increases with growing magnetic field, i.e., 
with growing n while that of the lateral velocity u(z), of the temperature 6(z), and of the 
magnetic field <p(z) is decreasing with increasing n without much change in the profiles. The 
discontinuous slope of 6{z) at z = ±| comes from the continuity of the heatflux. Here we 
used the normalization J dzw(z) 6(z) = 1. 

From the neutral stability problem one finds that the instability is a stationary one as 
for boundaries with infinite conductivity, ( = oo. The critical Rayleigh number r c (n) and 
the corresponding critical wave number k c (n) as function of the magnetic Rayleigh number 
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n are shown in Figs. [| and |], respectively, for 5 = 1 and several (. Both, k c as well as r c 
shift to smaller values when the conductivity of the plates decreases. The variation of r c (n) 
with n can be fitted reasonably well by a linear law. The fit coefficients are presented in 
Table |. for several values of (. 

The symbols in Figs. [3] and || refer to Rayleigh numbers and wave numbers that were 
observed at onset of convection in experiments done with fluid layers of different heights 
d betweeen 1 and 4 mm and an upper plate of fixed thickness dp = 6 mm with reduced 
conductivity £ = 5 ||. Our linear analysis showed for such values of S > 3 that r c and k c 
had well reached their large-5 asymptotic plateau behaviour: Both, r c and k c decrease in 
the range < 5 < 1 with growing 5 and have approached their large-5 constant asymptote 
around 5 = 1 so that at 5 > 1 there is no more significant variation with 5. Note in addition 
that the influence of S on convection properties decreases with growing conductivity ( of 
the plate — the thermal stress on the fluid, ATp, that determines the driving (|1.2|) and 
( |1.3| ) approaches AT for £ — > oo and the fluid approaches the limiting situation of beeing 
bounded by perfectly heat conducting plates. 

Our results agree in the nonmagnetic limit, n = 0, with the linear stability analysis of 
Jenkins et al. j3J for all (. Furthermore, in the presense of a magnetic field we confirm the 
theoretical results of Schwab [|J and Stiles et al. reported for the limit ( = oo. We have 
no explanation for the fact that the wave numbers (symbols in Fig. ^) that were observed 
to be selected by convective patterns close to onset in the experimental setup of Schwab et 
al. 0] lie closer to the ( = 1 theoretical critical wave numbers than to the ( = 5 curve. 

IV. WEAKLY NONLINEAR ANALYSIS 

Shortly above the stability threshold the nonlinear convective state can approximately 
be described by 

Y(x, z, t) = (u, w, 6, <p) T = A(x, t) Y(z) e ikcX + c.c. (4.1) 
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Here Y(z) is the eigenvector of the linear problem at the critical point presented in Fig. ^ and 
A(x, t) is the saturation amplitude of the critical mode determined by the nonlinearity in the 
field equations. For slightly supercritical driving (0 < e = r-r c C 1) a small band of lateral 
wave numbers with \k — k c \ = 0{e%) is excited giving rise to a slow spatiotemporal variation 
of the convective roll structure. These slow variations of the amplitude A are conveniently 
captured by the method of multiple scales |7| . It allows to derive in a straightforward manner 
the Ginzburg-Landau amplitude equation || 

T d t A(x,t)=[e + e d 2 x -l\A(x,t)\ 2 ]A(x,t) (4.2) 

from the hydrodynamic field equations. 

We have determined the coefficients r , £q, and 7 of this equation for rigid boundaries 
with different £ by numerically evaluating special scalar products of the eigenfunctions || . 
Since the coefficients with no magnetic field, i.e., for an ordinary fluid have not yet been 
presented in the literature as a function of ( we show them in Fig. |5] for 5=1. The linear 
coefficients r and £q are presented in Fig. |6| as functions of the magnetic Rayleigh number 
n. The nonlinear coefficient 7 shows a linear dependence on n, i.e., 

7(C,n) = 7o(C) + MCK (4-3) 

The expansion coefficients 70(C) an d & 7 (C) are given in Table [I] for several (. 

To evaluate the slope of the Nusselt number at the critical point we start from the 
relation 

Nu(r, n) 

that holds close to onset, r — r c {n) <^ 1 

dNu 
dn c 

Since r c (n) ~ 1 — b r n one can rewrite ([ 
more compact form 



1 + r _ZlM (4.4) 

7(n) 



Thus 

1 dr c (n) 
7(71) dn 



(4.5) 



as 7 = 70 + b 7 (l — r c )/b r to obtain ^f| c in the 



dNu 
dn 



1 

a — br r 



(4.6) 



The fit coefficients b = p- and a 



b + p are presented in Table |I[ 



V. PRANDTL NUMBER DEPENDENCE 

The results presented so far refer to fluids of Prandtl number a = 1. Now we discuss the 
a dependence for the case of perfectly heat conducting plates, ( = oo. The a dependence 
for finite (, say, £ > 1 is similar. The marginal stability curve and the linear eigenf unctions 
of the stationary eigenvalue problem have no Prandtl number dependence. Thus r c , k c , and 
£o are cx-independent. But r depends on a. Graphs for different n are given in Fig. |8] for 
( = oo. Note that the difference To(<7, n) — To(a = l,n) is practically independent of n. This 
implies that r (a,n) can be approximated well by 

cr + 0.5117 . , 

r (a,n) - r (a = l,n ) = — — 0.0769 5.1 

19.65 a 

for C = oo and the right hand side of ( |5.1|) is just r (cr, n = 0) [JTOj for the case without 
magnetic field. Thus one only needs to know t (o~ = l,n) which can be taken from the 
long-dashed curve in the lower part of Fig. [|. It represents ro(a = l,n,( = 1000) which is 
practically the same as To(a = 1, n, £ = oo). 

The a dependence of the nonlinear coefficient 7 can be fitted by the expression 

jM ^^ Hn) .fM + fM, (5.2) 

The fit coefficients 7^(71) are presented in Table |T| for several n. For a = 1 the n variation 
is practically linear, i.e., 7(72, o = 1) = 0.7027 + 1.42427 n. Graphs of the difference 7(71, a) — 
7(71, <T = 1) are given in Fig. || For Prandtl number larger than 10 the dependence on a can 
be neglected and one can use a = 00. 

Support by the Deutsche Forschungsgemeinschaft is gratefully acknowledged. 
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TABLES 



c 


0.8 


1 


2 


5 


10 


100 


1000 


oo 


r c = 1 — b r n 


Ra° c 


1256.9669 


1303.4412 


1439.8249 


1574.0344 


1634.8154 


1699.8132 


1706.9597 


1707.7618 


b r 


2.424 


2.403 


2.349 


2.304 


2.283 


2.265 


2.262 


2.259 


7 = 7o + b 1 n 


7o 


0.6120 


0.6145 


0.6330 


0.6634 


0.6802 


0.6991 


0.7023 


0.7027 




1.0709 


1.1223 


1.2656 


1.3883 


1.4468 


1.4992 


1.4594 


1.4243 


dNu 1 ._. 1 
dn Ic a—br c 


a 


0.434 


0.450 


0.498 


0.549 


0.575 


0.601 


0.596 


0.590 


b 


0.182 


0.194 


0.229 


0.262 


0.277 


0.292 


0.285 


0.279 



TABLE I. Fit-coefficients for the n dependence of r c ,7, and for different £. 



n 


7 (o) 


7 (2) 


7 (L 





0.699 


0.00832 


0.00471 


0.01 


0.710 


0.00966 


0.00453 


0.04 


0.745 


0.00941 


0.00547 


0.08 


0.792 


0.00916 


0.00673 


0.2 


0.945 


0.00851 


0.00102 


0.4 


1.212 


0.00728 


0.01411 


0.8 


1.756 


0.00568 


0.01867 


1 


2.014 


0.00502 


0.02035 



TABLE II. Fit coefficients for the a dependence of j(n,a) = 7<°)(n) - :L -f^ + for 
different n. 
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FIGURES 

FIG. 1. Vertical temperature profile of the quiescent conductive state in the plates (P) and 
in the fluid (F) here shown for 5 = ^ = 1. With AT beeing imposed externally the conductive 
temperature difference across the fluid layer is ATp = ±+J/£ ■ 

FIG. 2. Vertical variation of the linear critical eigenfunctions of lateral velocity u, vertical 
velocity w, magnetic field (p, and temperature 9 for ( = 5, 6 = 1, and different reduced magnetic 
Rayleigh numbers n. 

FIG. 3. Critical Rayleigh number r c versus magnetic Rayleigh number n for different £ and 
5 = l,x = 1- The curves for £ = 5 and £ = oo are practically identical. Symbols refer to 
experimentally observed convective thresholds [3] in setups with £ = 5 and 5 = 3(v), 4(o), 6(A), 
8(d), and 12(o). See text for further discussion. 

FIG. 4. Critical wave number k c versus magnetic Rayleigh number n for different £ and x = 1- 
Symbols refer to experimentally observed wave numbers of convection patterns close to threshold 
[3] in setups with ( = 5 and different 6 identified in Fig. ||[ See text for further discussion. 

FIG. 5. Coefficients of the amplitude equation as functions of £ without magnetic field, n = 0. 
The normalization of the critical eigenfunction entering into 7 was chosen to be J dzw(z) 6(z) = 1. 

FIG. 6. Linear coefficients tq and £q of the amplitude equation as functions of the magnetic 
Rayleigh number n for different £. The long-dashed curves for £ = 1000 are practically the same 
as for perfectly heat conducting plates with ( = 00. 

FIG. 7. Initial slope of the Nusselt number versus critical Rayleigh number. Symbols refer to 
experimental results obtained with plates of large heat conductivity [3]. 

FIG. 8. Linear coefficient To(a,n) — to(<t = 1, n) versus a for different magnetic Rayleigh 
numbers n. 
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FIG. 9. Nonlinear coefficient 7(0", n) — 7(0- = l,n) versus cr for different magnetic Rayleigh 
numbers n. 
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